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INFERRING  FINITE-TIME  PERFORMANCE 
IN  THE  M/G/l  QUEUE I NG  MODEL 

P . A .  Jacobs 
D.P.  Gaver 

1  .   INTRODUCTION  AND  SUMMARY 

The  M/G/l  model  for  a  single  server  approached  by  a  Poisson  (A)  stream 
of  arrivals  with  independent  but  otherwise  arbitrarily  distributed  service 
times  is  a  standard  tool  in  operations  research.    It  has  been  used  to 
model  situations  occurring  in  road  traffic,  computer  and  communication 
performance  evaluation  and  the  military,  and  has  been  modified  to 
accommodate  priority  and  server  breakdown  situations,  as  well  as  both 
heavy  and  light  traffic  situations  (p    =  AE  [S]  either  near,  but  below. 
unity,  and  near,  but  above ,  zero). 

Although  the  M/G/l  model  is  rather  simplistic,  little  beyond  the 
moments  of  its  long-run  or  steady-state  distribution  are  easily  available 
in  anything  like  closed  (and  simple)  analytical  form  readily  evaluated 
numerically'.    In  the  case  of  the  exponential  service  time  M/M/l  system, 
the  complete  analytical  solution  for  transition  probab i 1 i st i es  in  terms  of 
series  of  Bessel  functions  has  been  long  known,  and  numerical  transform 
inversion  by  Middleton  (1979)  has  provided  some  useful  tables  describing 
the  manner  in  which  steady  state  is  approached.   Recent  work  of  Abate  and 
Whitt  (1989,1988)  provides  some  handy  approximations,  particularly  in  the 
M/M/l  case.   The  work  of  Asmussen  and  Thorisson  (1988)  and  Newell  (1971, 


19S2)  furnishes  diffusion  approximations  and  heavy  traffic  results  as 
well. 

In  this  paper  we  explore  the  following  operationally  relevant 

Probl em :   A  single  server  is  approached  by  a  stream  of  arrivals  to  be 
modelled  as  Poisson.   One  has  available  a  sample  of  service  times  obtained 
by  observing  the  system.   Using  the  arrival  rate  (assumed  known,  at  least 
initially)  and  the  observed  service  times,  estimate  the  expected  waiting 
time  encountered  by  a  new  arriving  customer  at  finite  time  t.   Further- 
more, provide  an  assessment  of  estimation  uncertainty,  e.g.  confidence 
1 im  i t s  . 

The  above  is  just  one  prototypical  question  to  be  asked--it  is  the  one 
addressed  in  this  paper.   Other  questions  could  concern  the  probability 
that  an  arrival  at  t  will  wait  for  at  least  time  w,  or  that  the  customer 
waiting  time  has  never  exceeded  w  in  time  t  (starting  from  an  empty 
system) .   There  are  many  other  such  measures  of  system  performance  that 
are  relevant.   We  concentrate  here  on  inferring  the  mean  waiting  time  at 
finite  time  t,  and  do  so  non -paramet r  i  cal ly ,  i.e.  without  specifying,  and 
estimating  parameters  in  a  small  conventional  parametric  family  such  as 
the  Gamma.   Our  approach  is  to  work  with  the  well-known  Laplace  trans- 
form of  the  mean  (virtual)  waiting  time  for  the  M/G/l  system:   cf .  Takacs 
(1962)  and  Gaver  and  Jacobs  (1987),  an  empi  r ical  version  of  which  can  be 
obtained  by  utilizing  the  empirical  transform  of  the  busy  period,  itself 
depending  upon  the  empirical  transform  of  the  service  time  df  ,  derived 
from  the  observed  data. 

We  look  upon  our  present  results  as  exploratory  and  tentative,  but  of 
an  accuracy  useful  until  refinements  are  available. 


A  recent  paper  by  Tengels  (1988)  provides  mathematically  rigorous 
asymptotic  (large  sample  size)  treatment  of  certain  aspects  of  our 
nonparamet r i c  inference  problem.   Application  of  bootstrap  methods,  cf . 
Efron  (1979)  can  yield  further  insights  into  the  inference  at  the  price  of 
computer  intensivity. 


2.   MATHEMATICAL  FORMULATION 

Consider  an  M/G/l  queue  with  known  arrival  rate  A  and  independent 
identically  distributed  service  times,  S,  having  nth  moment  E  [Sn]  <  oc  . 
Suppose  there  are  no  customers  either  waiting  or  being  served  at  time  0. 

Let  to't  be  the  virtual  waiting  time  at  time  t  and  put  (conditional  on 
W0  =  0) 


VW  (  s  )  = 


e~st"E[Wt]  dt,    s  >  0 


(1) 


We  know  that 


/%     p    -    1  /    \     1 

vw(s)  =  —    +  Poo(s)  s  ; 


(2) 


where : 


PooC5;  = 


s  +  A[l  -  bi^s)]  ' 


(3) 


b(s)  being  the  Laplace  transform  of  a  busy  period  duration;  b(s)  is  the 
smallest  positive  solution  of 

b(s)  =  Fs(s  +  A(l  -  b(s)))  ; 


(4) 


Fs  is  the  Lapl ace-St i e 1 t j es  transform  of  the  service  time.   See  Takacs 
(1962)  and  Gaver  and  Jacobs  (1986) . 

Two  approximations  to  E [W»  |  W«  =  0]  will  be  detailed  below;  one  for  the 
case  of  the  traffic  intensity  AE  [S]  <  1;  the  other  for  the  case  of  the 
traffic  intensity  AE  [S]  >  1.   Both  are  based  on  the  transforms  (2),  (3), 
and  (4)  . 


3.   APPROXIMATIONS  FOR  STABLE  QUEUES:   EXPONENTIAL  APPROACH 

In  this  subsection,   it  is  assumed  that  p    -    AE  [S]  <  1.    In  this  cas< 

there  is  a  limiting  expected  virtual  waiting  time 

:2- 


C  i  •  rrwl  AE[S^1 

U-    =    \l?x    E^tl     =    2(1    -    p)     ' 


(5) 


with 


p    =    AE[S]  , 


(6) 


the  traffic  intensity. 

Our  choice  of  the  form  of  the  approximation  is  based  on  an  observation 
of  Odoni  and  Roth  (19S3) ,  and  on  relaxation  time  calculations  made  by 
Keilson  (1979),  Newell  (1971)  and  others.   The  explicit  approximation  of 
the  exponential : 

E[Wt]  *  £00(1  -  e_it) 


X1EfWt]  %  e-^t 


(7) 


1  oc 


It  should  be  noted  that  in  many  cases  the  exponential  approach  to  the 

limiting  value  is  not  exact,  being  only  an  approximation,  cf .  Asmussen  and 
Thorissan  (1988) .   Setting 


se  st    u  t-_ 2£  dt  = 


0 


-st   -^t  _> 
se  Sl  e    dt 


(S) 


we  obtain 


STi>w(s)  -  W- 


Uoc  s  +  0 


(9) 


which  is  now  solved  for  (3    at  selected  s-values: 


i  +  Hh)  = 


svw(s) 


nb  = 


F  i  nal 1 v 


0    = 


•woc    -     [g^w(s)    -    ^3c] 

sVw  (s)     ~     "oc 


^w(s) 


(10) 


st-w  (s)     -    Woo 
The    approximation    to    E[Vt]     is 

E[Wt]  *  WooCl  -  e'0t) 


(11) 


(12) 


where  0  is  evaluated  as  above  at  s  =  -.  It  turns  out  that  the  ^-values  so 
obtained  are  not  constant  but  are  (extremely)  slowly-changing  functions  of 
t.  as  should  be  the  case  from  asymptotic  analysis. 

Table  1  reports  the  values  of  the  approximating  E  [W't]  for  an  M/M/l 
queue  with  A  =  1  and  E  [S]  =  .95.   The  approximating  values  are  compared  to 
values  obtained  by  Middleton  (1979) ,  who  used  a  numerical  inversion  of  the 
Laplace  transform  V'w  due  to  Stehfest  (1970) .   The  approximation  is  very 
good  for  the  times  less  than  or  equal  to  100.   For  the  times  greater  than 
100  the  approximation  is  within  11%  of  the  true,  and  conservatively 
slightly  overestimate  of  the  true  mean. 


Table  1 

Exponential  Approximation 
M/M/l  Queue:   A  =  1,  E [S]  =  .95 


Time 

True 

t 

E[Wt] 

20 

3.9 

40 

5.5 

60 

6.6 

80 

7.4 

100 

8.1 

400 

13.2 

800 

15.5 

1200 

16.6 

Note  : 

\>dc 

=    18.05 

Exponent  i  al 

Approxi mat  ion 

3.7 

5.4 

6.6 

7.6 

8.4 
14.6 
17.0 
17.  7 


Table  2  reports  the  values  of  the  approximating  E  [Vt]  "for  an  M/G/l 
queue  with  A  =  1  and  gamma  service  times  with  shape  parameter  q  =  .2  and 
scale  parameter  3    =    4.5.   Once  again  the  approximation  is  very  good.   The 
values  of  the  true  E  [Wt]  are  from  Middleton  (1979)  . 


Table  2 

Exponential  Approximation 

M/G/l  Queue:   A  =  1 

Service  Time  Distribution  Gamma  shape  parameter 

q  =  .2.  scale  parameter  0    =  4.5,  E [S]  =  .90 


Time  True                    Exponential 

t  E  [Wt]                   Approximation 

4.9  2.3                        2.2 

9.7  3.7                        3.4 

29.2  6.9                        6.7 

48.6  8.8                                                           8.8 

97.2  11.9  12.5 

194.4  15.3  16.9 

486  20.1  21.9 

972  22.6  23.7 

1458  23.5  23.8 

1944  23.9  24.1 

Note:   Woe  =  24.3 


The  numerical  results  obtained  suggest  that  in  general  the 
approximation  is  biased  slightly  on  the  high  side  for  large  t.   It  is  thus 
conservative  in  the  sense  that  predictions  of  expected  waiting  times  tend 
to  be  slightly  overstated. 
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4.   THE  HEAVY  TRAFFIC  EXPONENTIAL  APPROXIMATION  (HTE) 

The  most  computational ly-intensive  part  of  the  exponential 
approximation  is  the  numerical  evaluation  of  the  transform  of  the  busy 
period  b(s)  using  equation  (4).   This  computation  has  been  done  by  search. 
In  this  section  another  approximation  is  described  which  uses  a  heavy 
traffic  approximation  to  obtain  b(s),  and  is  computationally  less  arduous. 

Let  {B(t)}  be  a  Brownian  motion  with  drift  v    =  AE  [S]  -  1  and 
infinitesimal  variance  a      =  AE  [S  ]  where  S  is  a  generic  service  time.   Let 
Tx  be  the  first  passage  time  to  state  0  given  B(0)  =  x.   Put 


*(s;x)  =  E[e"sTx] 


(13) 


It  is  known  that 


© (s : x)  =  e 


q(s)x 


(14) 


where 


a(s)  = 


|^2  +  2o-2s 


(15) 


Of  course  this  function  is  explicit  and  readily  evaluated,  and  no  search 
i  s  needed  . 

Approximating  the  virtual  waiting  time  process  by  the  Brownian  motion 
with  drift,  the  proposed  approximation  for  the  transform  of  the  busy- 
period  is 


bBM(s)  =  E[ea(s)S]   , 


(16) 


the  Laplace-St ie 1 j es  transform  of  the  service  time  distribution  evaluated 
at  o(s) .   The  remainder  of  the  approximation  is  the  same  as  the 
exponential  approximation.   Table  3  shows  the  values  of  the  heavy  traffic 


approximation  for  the  M/G/l  queue  with  A  =  1  and  gamma  service  time 
distribution  with  shape  parameter  a  =  0.2  (very  long-tailed)  and  scale 
parameter  ,3    =    4.5.   The  heavy  traffic  approximation  is  not  as  good  as  the 
exponential  one  and  tends  to  overestimate  E  [Wt]  .   However,  the 
approximation  is  still  practically  adequate,  and  is  far  easier  to  compute 
than  is  the  "exact"  exponential  approximation  of  Section  2. 

Table  3 

Unsaturated  Heavy  Traffic  Approximation 

M/G/l  Queue:   A  =  1,  Gamma  Service  Times 

Shape  parameter  a  =  0.2,  scale  parameter  0    =  4.5,  E [S]  =  .90 

Heavy  Traffic 
Approx  imat  ion 

2.5 

3.9 

7  .  5 

9.8 

13.6 

17.9 

24.  1 

24.3 


T  ime 

True 

t 

E  [Wt] 

4.9 

2.3 

9.7 

3.7 

29.2 

6.9 

4S.6 

8.8 

97.2 

11  .9 

194.4 

15.3 

972 

22.6 

L944 

23.9 

s'ote  : 

Vac 

=    24 

.3 
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5.   THE  BUSY  PERIOD  APPROXIMATION 

This  approach  is  based  on  approximating  the  rate  of  approach  of  the 
probability  o  an  empty  system,  PqoC^)'  *-°     its  eventual  value.   Since  the 
process  of  alternating  busy  and  idel  periods  can  be  expressed  in  terms  of 
an  alternating  renewal  process,  methods  of  Gaver  and  Jacobs  (1988)  can  be 
appl i  ed . 


Note 


s  •  tfw(s)  =  P— g —  +  p00(s)  ,  (17) 


where 


Poo(s)  =  s  j1.  b,s  (IS) 

1   f   AL S J 

and 

Poo(s)  =      e"st  P{\(t)  =  0|X(0)  =  0}dt  (19) 

where  X(t)  is  the  number  of  customers  waiting  or  being  served  at  time  t. 
Assume  the  queue  is  stable;  then 

lirn^  P{X(t)  =  0|X(0)  =  0}  =  E  [Id  ^f  ^[Busy]  =  1  "  >  '  <20> 

where  Idle  is  the  length  of  an  idle  period  and  Busy  is  the  length  of  a 

busy  period.   Approximate  as  follows  focussing  on  the  probability-  that  the 
system  is  empty,  approximates  as  follows  in  terms  of  the  parameter  3b  : 

P{X(t)  =  0|X(0)  =  0}  «  (1  -  p)    +    pe~0t>t     ;  (21) 

from    this    comes 

E[Wt]     «     (p    -    l)t    +     (1     -    p)t    +    -f  (1     -    e^bt) 

=    f  (1     -    e"V)     .  (22) 

11 


To  find  /?h: 


Poo(s)  =  * s-^  +  0     \    s  (23) 


poo(s)  -  I  ~    ,,11-^J.  ■  i  l  +^c(s)  (24) 


1  +  AC- s11"^] 


where 


c<*>  -  [1  ;E:bB?)]     "d     ^m  =  H-y  <25; 

Solving  for  0b    results  in 


sC(s)[r4-7  +  1]    sC(s)  [^] 
1  -  C(s)      ~    1  -  C(s) 


3b    ~    i F7^ -  i — — F7T^ '  (26) 


and  evaluate  at  s  =  ;* 


,  1 
t  ' 

Table  4  shows  values  for  the  busy  period  approximation  for  an  M/M/l 
queue  with  A  =  1  and  p    =  .95.   The  approximate  values  are  below  the  true 

values  for  most  times  t.   This  approximation  is  not  as  accurate  as  is  the 
exponential  approximation  of  Section  2,  no  longer  being  conservative, 

Table  4 

Busy  Period  Approximation 

M/M/l  Queue 

A  =  1  ,  E[S]  =  .95 

Time                   True  Approximate 

t                   E  [Wt]  E  [Wt] 

20                    3.9  4.1 

40                    5.5  5.4 

60                    6.6  6.3 

80                    7.4  7.0 

100                    8.1  7.5 

400                   13.2  11.6 

800                   15.5  13.5 

1200                   16.6  14.6 
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6.   THE  INSTABLE  OR  SATURATED  QUEUE  EXPONENTIAL  APPROXIMATION 
Rewriting,  we  have 


=  L^l   +  i 


s2  -  -  s  s  +  AQl  -  b(s)] 


(25; 


In  this  subsection  we  assume  that  p    >  1.   When  this  is  true  the 
distribution  of  the  length  of  the  busy  period  is  not  honest,  and 
b(0)  <  1.   We  will  write  b(0)  b#(s)  =  b(0)  [grrrh  ,  where  b#(s)  is  now  the 
transform  of  an  honest  random  variable.   Now 

vw(s)  =  P—I-1     +  1   1 .  ,'29; 

s^        s  +  A[l-  b(0)  b*(s)] 

For  smal 1  s 

♦"(•)  *  '-ir  *  l  *[i  -Vo.j  ■  (30) 

Thus  for  large  t 

E[Wt]  *  (p  -  l)t  +  A[1  _lb(0)]  .  (31) 


Furthe  r 


Vw(s)  =  ^^  +  |o(s)  (32) 


where 


(s)  = 


s  +  A[l  -  b(0)b#(s)]    s  +  A  -  Ab(0)b#(s) 
Multiplying  through  by  s  +  A  yields 

(s  +  A)<Ks)  = 


i  -  ;rf-iMO)b*(s) 


or 
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(s  +  A)o(s)  =  1  +  (s  +  A)e(s)|"s  *  Ab(0)b#(s)j 


(33) 


and 


o(s)  =  j-L-j  +  y-J-j  b(0)b*(s)*(s) 


(34) 


rOC 


Lett  i  ng  <p  (s)  = 


e'st  f(t)dt 


f(t)  =  e"At  + 


f(t-u)h(u)du 


where 


7(s)  = 


oo 


e"sth(t)dr  =  b(0)  j-A-j  b#(s) 


(35) 


and 


h(0)  =  b(0)  <  1  . 


Thus  f    satisfies  a  renewal-type  equation  with  a  degenerate  distribution 
and  it  can  be  shown  ( c"f  .  Feller  (1966)) 

-0r. 


f(t)  *  ce 


as     t  —  oc 


(36) 


Thus  .  since 


^W ( s )  - 


P-^    +  I  *(■) 


E[Wt]  *  (p    -    l)t  +  S(l  -  e_/?t) 


(37) 


as  t  — •  oo .   Comparing  (36)  with  (31)  indicates  that 


c  _ 


p  A[l  -  b(0)] 


(38) 
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In  the  saturated  case  the  exponential  approximation  is  computed  a« 
foil ows : 

E[Wt]  -  (p    -     l)t  *  a  (1  -  e"5t) 


where 


a  = 


A[l  -  b(0)] 


Thus 


E[Wt]  -  (p    -  l)t  -  a      -0t 


(39) 


Taking  Laplace  transforms 


s^w(s)  ~  (— s — )  "  a 


i  +  /?(i) 


(40) 


which  leads  to 


0    = 


-sz^w(s)  +  (p  -  l) 
svw(s)   -   (£-^ )   -  1 


(41) 


s  i  nee 


a  = 


\  [  i  _  b  (  0  )  ] 


(42) 
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7.   THE  SATURATED  QUEUE  HEAVY  TRAFFIC  EXPONENTIAL  APPROXIMATION 

Once  again  the  most  arduous  computation  in  the  saturated  queue 
exponential  approximation  is  the  computation  of  b(s) ,  the  transform  of  the 
length  of  the  busy  period.   As  before,  an  approximation  is  proposed  which 
uses  a  heavy  traffic  approximation  to  evaluate  b(s)  and  b(0) . 
Specifically,  since  v    —    AE [S]  -  1  >  0  let 

Q(0)  =  z2»  -  (43) 

a 

and 


a(s)  =  ^ — _ (44) 


a 


and  put 


and 


b(0)  =  E[eQ(0)s]  (45) 


b(s)  =  E[eQ(s)s]   .  (46) 


The  remainder  of  the  approximation  is  as  in  the  exponential  approximation 
for  the  saturated  queue. 

Table  5  shows  results  for  both  the  "exact"  exponential  approximation 
and  heavy  traffic  approximation  for  a  saturated  M/G/l  queue  with  A  =  1  and 
gamma  service  times  having  shape  parameter  0.2  and  scale  parameter  7.5  and 
E [S]  =  1.5.   The  true  values  are  from  Middleton  (1979).   As  is  expected 
the  exponential  approximation  gives  values  which  are  closer  to  the  true. 
Both  approximations  improve  as  t  increases.   The  exponential  approximation 
is  very  good  for  all  times  except  the  very  smallest. 
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Table  5 

Saturated  M/G/l  Queue 

A  =  1.  Gamma  Service  Times,  E [S]  =  1.5 

Shape  parameter  0.2,  Scale  parameter  7.5 


Time 

True 

Exponent  ial 

Heavy    Traffic 

t 

E[Wt] 

App 

roximat  i  on 

App r ox  imat  ion 

1  .6 

3.9 

1  .9 

2.  1 

8.1 

7.5 

7.4 

8.2 

16.2 

12.9 

13.  1 

14.4 

27.0 

19.3 

19.7 

21  .6 

43.2 

28.3 

28.8 

31  .3 

54.0 

34.0 

34.6 

37.3 

108.0 

61  .9 

62.2 

65.6 

162.0 

89.  1 

89.2 

92.8 

216.0 

116.2 

116.2 

119.8 

324.0 

170.2 

170.2 

173.8 

432.0 

224.2 

224.2 

227.8 

8.   INFERENCE 

The  approximations  to  E[Wt]  can  be  used  for  inferential  purposes  by 
replacing  moment  and  transforms  by  sample  moments  and  empirical 
transforms.   More  specifically,  suppose  A  =  1  is  known  and  service  time 
data  dj^  ,  .  .  .  ,dn  are  collected.   The  empirical  Laplace  transform  of  the 
service  time  distribution  is 

n    -sd: 
*s(s)  =  E  e   '  ;•  (47) 

i=l 

^'oc  is  estimated  by 


where 


d  =  n  £  d,  ;  (49) 

i=l 

d5  =  h    t     d,2  ;  (50) 

i=l 


and 


P    =  Ad  .  (51) 

The  empirical  Laplace  transform  of  the  service  times  and  the  moment 
estimator  of  p,     p,     can  be  used  in  formulas  (2)-(4)  to  obtain  an  estimate 
of  vw(s)  •   This  estimate,  Vw(s)  »  "together  with  £  can  be  used  in  the 
approximations  to  obtain  estimates  of  E[W(t)] . 

Tables  6-8  report  results  of  simulation  experiments  to  study  the 
behavior  of  the  estimates.   Each  simulation  has  1000  replications.   Each 
replication  supposes  a  sample  of  100  service  times.   The  random  numbers 
were  generated  using  LLRAND0MII  random  number  package,  see  Lewis  and  Uribe 
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Table  6 

Moments  of  Estimates  of  E[Wt' 

for 

M/M/l  Queue  with 

A  =  1   E[S]  =  0.9 


Time 

t 

6, 

.5 

32 

.4 

97, 

.2 

162. 

0 

486, 

,0 

64S, 

,0 

True 

E[Wt] 

1  .9 
4.2 
5.9 
6.7 
7.9 
8.0 


Exponent  i  al 

\pproxi mat  ion 

Mean 

Var  i  ance 

1  .8 

0.  1 

4.3 

2.0 

6.8 

14.0 

8.3 

33.1 

13.0 

210.6 

14.7 

347.6 

Heavy 

Traf f  i  c 

Approxi mat  ion 

Mean 

Var i  ance 

2.0 

0.  1 

4.5 

2.4 

7.  1 

15.6 

8.6 

35.9 

13.2 

218.5 

14.9 

357  .  5 

Table  7 

Moments  of  Estimates  of  E [W't] 

M/G/l  Queue 

A  =  1   E[S]  =  0.5 

Gamma  Service  Times 

Shape  =  0.2    Scale  =  2.5 


T  ime 

t 

0.3 
1  .8 
3.0 
4.8 
9.0 
24.0 


True 

E[Wt] 

0.13 
0.53 
0.71 
0.90 
1  .2 
1  .4 


Exponent  i  al 

Approx  imat  ion 

Mean 

Variance 

0.  12 

0.001 

0.51 

0.024 

0.70 

0.058 

0.91 

0.13 

1  .2 

0.33 

1  .6 

1  .0 

Heavy 

T 

raff  i  c 

Approx  imat  ion 

Mean 

V 

ar i  ance 

0.13 

.001 

0.56 

.03 

0.78 

.08 

1  .0 

.  17 

1  .3 

.44 

1  .6 

1 

.2 
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Table  8 

Moments  of  Estimates  of  E [W.] 

M/G/l  Queue 

A  =  1  ,   E[S]  =  1  .5 

Gamma  Service  Times 

Shape  =  0.2    Scale  =  7.5 


Time 

True 

Exponent  ial 

Heavy 

Traf f  i  c 

t 

E[Wt] 

Approximat  ion 

Approximat  ion 

Mean 

Variance 

Mean 

Vari ance 

1  .6 

3.9 

1  .9 

0.30 

2.0 

.31 

8.1 

7.5 

7.3 

6.7 

8.1 

7.6 

54.0 

34.0 

34.5 

256.8 

36.8 

277.7 

162.0 

89.1 

90.1 

2436.0 

92.9 

2488.0 

324.0 

170.2 

171  .3 

10,210 

174.3 

10.296 

(1981).   If  the  true  average  service  time  is  close  to  1,  then  it  is 
clearly  possible  for  the  sample  traffic  intensity  p    to  be  less  than  or 
greater  than  1.    If  p     is  less  than  1,  then  the  estimate  using  the 
approximation  for  the  stable  queue  is  computed.   If  p    >  1,  then  the 
estimate  using  the  approximation  for  the  unstable  queue  is  computed.   This 
choice  appears  to  be  natural  unless  other  information  is  available,  or 
more  assumptions  are  made. 

In  Tables  6-8  are  reported  means  and  variances  of  the  estimates  of 
E[Wt]  for  the  estimates  based  on  the  exponential  approximation  and  the 
heavy  traffic  approximation. 

In  Table  6  results  for  a  M/M/l  queue  with  E [S]  =  .90  are  given.   The 
same  random  numbers  were  used  to  compute  the  estimates  for  each  of  the 
times  t.   Of  the  1000  replications,  137  had  p    >    1  and  so  the  unstable 
queueing  approximations  were  used  in  these  cases.   Doubtless  it  is  the 
contribution  of  these  cases  that  lead  to  the  pronounced  over  estimate  of 
the  mean . 
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In  Table  7  are  reported  results  for  an  M/G/l  queue  with  A  =  1  and 
gamma  service  times  with  E  [S]  =  .5  having  shape  parameter  .2  and  scale 
parameter  2.5.   None  of  the  1000  replications  has  p    >  1. 

In  Table  S  are  reported  moments  for  the  estimates  of  E[Vt]  for  ar. 
unstable  M/G/l  queue  with  A  =  1  and  gamma  service  times  with  shape 
parameter  0.2  and  scale  parameter  7.5  and  E [S]  =  1.5.   Of  the  1000 
replications  p    >  1  for  938.  of  them.   The  true  values  of  E[Wt]  in  all  the 
tables  are  from  Middleton  (1979). 

The  means  of  the  heavy  traffic  approximation  are  larger  than  those  for 
the  exponential  approximation.   The  variances  increase  as  t  increases. 
The  very  large  variances  in  Tables  6  and  8  are  attributable  to  those 
estimates  of  E [Wt]  for  which  p    >  1. 

In  order  to  put  the  large  variances  of  Table  8  into  perspective, 
consider  the  following  approximation. 

As  a  first  approximation,  if  AE[S]  >  1.  then  E [Wt]  %  (AE[S]  -  l)t 
which  can  be  estimated  using  the  mean  of  the  service  times 

m(t)  =  (Ad  -  l)t  . 
The  variance  of  m(t)  for  a  sample  size  of  100  is 
\ar[m(t)j  =  az        1QQ    **     . 

In  the  case  of  a  gamma  service  time  distribution  with  shape  parameter  0.2 
and  scale  parameter  7.5,  A  =  1,  t  =  162 

Var[m(t)]  =  2952.5  . 

Comparing  this  number  to  the  corresponding  variance  of  the  estimate  in 
Table  8  of  2436  indicates  that  the  latter  variance  is  not  unreasonable. 
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The  means  of  the  estimates  in  Table  7  are  close  to  the  theoretical 
values  in  this  stable  queue.   The  means  of  the  estimates  in  Table  6  are 
close  to  the  theoretical  values  for  the  smaller  t;s.   For  larger  t  the 
means  are  greater  than  the  theoretical  values  owing  to  those  replications 
for  wh  ich  p    >  1  . 

The  means  of  the  estimates  in  the  unstable  queue  case  of  Table  8  are 
close  to  the  true  values  of  E[Wt]  except  for  the  smallest  time,  t  =  1.6. 
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9.   CONCLUSIONS 

This  paper  proposes  an  easily-computed  approximation  to  the  finite- 
time  expected  waiting  time  "for  an  M/G/l  system  starting  from  an  empty 
condition.   Both  unsaturated  (p    <    1)  and  saturated  (p  >  1)  conditions  are 
considered.   Numerical  evidence  is  presented  to  indicate  that  the  quality 
of  the  approximation  is  usefully  good,  especially  when  ease  of  computation 
is  an  issue.   Further,  the  methodology  is  adapted  to  assess  expected 
waiting  time  when  inferences  must  be  made  from  a  random  sample  of  service 
times,  and  the  decision  is  made  to  do  so  nonparamet r i cal ly ,  i.e.  without 
fitting  a  specific  function.   The  results  appear  reasonable  and 
potentially  useful,  and  are  not  burdensome  to  obtain.   The  methodology 
investigated  can  also  be  applied  to  the  variety  of  queueing  models  that 
are  close  siblings  of  M/G/l:   priority  and  breakdowns  and  "vacations" 
being  examples.   Of  course  other  approximating  and  inferential  options 
remain  to  be  investigated. 


23 


References 

J.  Abate  and  W.  Whitt.   Transient  behavior  of  the  M/M/l  queue  via  Laplace 
transforms,  Adv .  Appl ■  Prob .  20  (1988),  pp.  145-178. 

J.  Abate  and  V.  Whitt.   Approximations  for  the  M/M/l  busy-period  distribu- 
tion.  Unpublished  paper,  1987. 

S.  Asmussen  and  H.  Thorisson.   Large  deviation  results  for  time-dependent 
queue  length  distributions.   Commun .  Statist. — Stochastic  Models  4 
(1988),  pp.  99-116. 

B.  Efron.   Bootstrap  methods:   another  look  at  the  jackknife.   Annals  of 
Stat  i  st  i  cs .  Vol.  7.  pp.  1-26. 

to" .  Feller.  An  Introduction  to  Probability  Theory  and  its  Applications. 
Vol.  II.   John  Wiley  and  Sons,  Inc.,  New  York,  1966. 

D.P.  Gaver  and  P. A.  Jacobs.   On  inference  and  transient  response  for  M/G/l 
models.   Teletraff ic  Analysis  and  Computer  Performance  Evaluation 
(Ed.  O.J.  Boxma.  J.W.  Cohen,  and  H.C.  Tijms),  Elsevier  Science 
Publishers  B.V.  (North-Holland),  1986. 

D.P.  Gaver  and  P. A.  Jacobs.   System  availability:   time  dependence  and 

statistical  inference  by  (semi)  non-parametric  methods.   Naval  Post- 
graduate School  Technical  Report,  1988. 

D.P.  Gaver.   Observing  stochastic  processes  and  approximate  transform 
inversion.   Operations  Research  15  (1966),  pp.  444-459. 

J.  Keilson.   Markov  Chain  Mode  1 s--Rar i ty  and  Exponent i al ity ,  Springer- 
Verlag,  New  York,  1979. 

P.A.W.  Lewis  and  L.  Uribe.   The  new  Naval  Postgraduate  School  random 

number  package --LLRAND0M  II.   Naval  Postgraduate  School  Technical 
Report  NPS55-81-005,  Monterey,  1981. 


24 


M.R.  Middleton.   Transient  effects  in  M/G/l  queues:   an  empirical  investi 

gation.   Technical  Report  No.  85,  Department  of  Operations  Research, 

Stanford  University,  Stanford,  Calif.  1979. 
A.R.  Odoni  and  E.  Roth.   An  empirical  investigation  of  the  transient 

behavior  of  stationary  queueing  systems.   Operations  Research  31 

(1983),  pp.  432-455. 
H.  Stehf est  ,  Algorithm  368-   Numerical  inversion  of  Laplace  transforms. 

Comm.  Assoc.  Comput .  Mach.  13  (1970),  pp.  47-49  [erratum  13,  624]. 
L.  Takacs.   Introduction  to  the  Theory  of  Queues.   Oxford  University 

Press,  New  York,  1962. 
G.F.  Newell.   Applications  of  Queueing  Theory.   Chapman  and  Hall,  London. 

1982  (Second  Edition). 


25 


DISTRIBUTION  LIST 

No.  of  Copies 

Library  (Code  0142)  2 

Naval  Postgraduate  School 
Monterey,  CA  93943-5000 

Defense  Technical  Information  Center  2 

Cameron  Station 
Alexandria,  VA  22314 

Office  of  Research  Administration  (Code  012)  1 

Naval  Postgraduate  School 
Monterey,  CA  93943-5000 

Center  for  Naval  Analyses  1 

4401  Ford  Avenue 
Alexandria,  VA  22302-C268 

Library  (Code  55)  1 

Naval  Postgraduate  School 
Monterey,  CA  93943-5000 

Operations  Research  Center,  Rm.  E40-164  1 

Massachusetts  Institute  of  Technology 
Attn:  R.  C.  Larson  and  J.  F.  Shapiro 
Cambridge,  MA  02139 

Koh  Peng  Kong  1 

0A  Branch,  DS0 
Ministry  of  Defense 
Blk  29  Middlesex  Road 
SINGAPORE  1024 

Arthur  P.  Hurter,  Jr.  1 

Professor  and  Chairman 

Dept  of  Industrial  Engineering 

and  Management  Sciences 
Northwestern  University 
Evanston,  IL  60201-9990 

Institute  for  Defense  Analysis  1 

1800  North  Beauregard 
Alexandria,  VA  22311 

Professor  H.  G.  Daellenbach  1 

Department  of  Operations  Research 
University  of  Canterbury 
Christchurch,  NEW  ZEALAND 

Donald  P.  Gaver  (Code  55Gv)  10 

Naval  Postgraduate  School 
Monterey,  CA  93943-5000 

Patricia  A.  Jacobs  (Code  55Jc)  10 

Naval  Postgraduate  School 
Monterey,  CA  93943-5000 


DUDUEYKNOXUBBARY 


§  2768  00337197  2 


